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Abstract. An experimental analysis of the asynchronous version of the 
"Game of Life" is performed to estimate how topology perturbations 
modify its evolution. We focus on the study of a phase transition from an 
"inactive-sparse phase" to a "labyrinth phase" and produce experimental 
data to quantify these changes as a function of the density of the initial 
configuration, the value of the synchrony rate, and the topology missing- 
link rate. An interpretation of the experimental results is given using the 
hypothesis that initial "germs" colonize the whole lattice and the validity 
of this hypothesis is tested. 



1 Introduction 

Cellular automata were originally introduced by von Neumann in order to study 
the logical properties of self- reproducing machines. Following Ulam's suggestions, 
the requirements he made for constructing such a machine was the discreteness 
of space using cells, discreteness of time using an external clock, the symmetry 
of the rules governing cells interaction, and the locality of these interactions; it 
resulted in the birth of the cellular automaton (CA) model. In order to make 
the self-reproduction not trivial he also required that the self-reproducing ma- 
chine should be computation-universal (e.g., [13]). The resulting CA used 29 
elementary states for each cell and updates used 5 neighbors. Later on, Conway 
introduced a CA called "Game of Life" or simply Life which was also proved to 
be computation universal [2]. This CA is simpler than von Neumann's in at least 
two ways: the local rule uses only two states and it can be summarized with by 
sub-rules (birth and death rules). 

However, the question remained open to know what is the importance of 
perfect synchrony on a CA behavior. Indeed, since the first study on the effects 
of asynchronous update carried out by Ingerson and Buvel [6], many criticisms 
have been addressed to the use of the CA as models of natural phenomena. 
Some authors investigated, using various techniques, how synchrony variations 
changed CA qualitative behavior [6], [10], [3], [5], [14], [8]. All studies agree on the 



fact that for some CA, there are situations in which small changes in the update 
method lead to qualitative changes of the evolution of the CA, thus showing the 
need for further studies of robustness to asynchronism. Similarly, some authors 
investigated the effect of perturbing the topology (i.e., the links between cells) 
in one dimension [11] by adding links, or in two dimensions with small-world 
construction algorithms [15], [9]. Here too, the studies showed that robustness 
to topology changes was a key factor in the CA theory and that some CA showed 
"phase transitions" when varying the intensity of the topology perturbation. 

The aim of this work is to question, in the case of Life, the importance of 
the two hypotheses used in the classical CA paradigm: what happens when the 
CA is no longer perfectly synchronous and when the topology is perturbed? In 
Section 2, we present the model and describe the qualitative behavior induced 
by the introduction of asynchronism and/or topology perturbations. In Section 
3.1, we observe that (i) Life is sensitive to asynchronism; (ii) robust to topology 
perturbations and (iii) that the robustness to asynchronism is increased when the 
topology characteristics become irregular. Section 3.2 is devoted to presenting a 
rigorous experimental validation and exploration of these phenomena for which 
a potential explanation based on the notion of "germ" development is discussed 
and studied in Section 4. 

2 The Model 

Classically, Life is run on a regular subset of 1? . For simulation purposes, the 
configurations are finite squares with N x N cells and the neighborhood of each 
cell is constituted of the cell itself and the 8 nearest neighbors (Moore neighbor- 
hood). We use periodic boundary conditions meaning that all cell position indices 
are taken in Z/NZ. The type of boundary conditions does play an important 
role at least for small configurations as shown in [4]. 

Life belongs to the outer-totalistic (e.g., [12], [11]) class of CA: the local 
transition rule / is specified as a function of the present state q(c) and of the 
number Si(c) of cells with states 1 in the neighborhood. The Life transition 
function f{q, S\) can be written: 

/(0, Si) = 1 if Si = 3 ; /(0, Si) = otherwise, (birth rule) 

/(l, Si) = 1 if Si = 2 or Si = 3;/(l, Si) = otherwise. (death rule) 

In the sequel, we consider Life as an asynchronous cellular automaton (ACA) 
acting on a possibly perturbed topology. 

There are several asynchronous dynamics: one may, for example update cells 
one by one in a fixed order from left to right and from bottom to top. This update 
method is called "line-by-line sweep" [14] and it has been shown that this type 
of dynamics introduce spurious behaviors due to the correlation between the 
spatial arrangement of the cells and the spatial ordering of the updates. These 
correlations can only be suppressed with a random updating of the cells. In 
this work, we choose to examine only one type of asynchronism which consists 



in applying the local rule, for each cell independently, with probability a. The 
parameter a is called the "synchronicity" [5] or the synchrony rate; one can also 
view it as a parameter that would control the evolution of a probabilistic cellular 
automata (PCA) where the transition function results in applying Life rule with 
a probability a and the identity rule with probability 1 — a. 

We choose to perturb topology by definitely removing links between cells. Let 
Go = (£, E ) be the oriented graph that represents cells interactions: (c, d) e 
E if and only if c' belongs to neighborhood of c. The graph with perturbed 
topology G = (C, E) is obtained by examining each cell c e £ and, for each 
cell in the neighborhood of c and removing the link (c, c') with a probability e~; 
the parameter e~ is called missing-link rate. Note that, as the local function is 
expressed in an outer-totalistic mode, we can still apply it on neighborhoods of 
various sizes. The definition we use induces an implicit choice of behavior in the 
case where a link is missing : the use of Si in the local rule definition implies 
that the cell will consider missing cells of the neighborhood as being in state 0. 
Other choices would have been possible; for example assuming this state to be 
1 or the current value of the cell itself. 

3 Observations and Measures 
3.1 Qualitative Observations 

Figure 1 shows that the behavior of Life depends on the synchrony rate a: a phase 
with labyrinthine shapes appears when a is lowered. Bersini and Detours studied 
this phenomenon and noticed that the asynchronous (sequential) updating of 
Life was significantly different from the (classical) synchronous version in that 
sense that a "labyrinth phase" (denoted by LP) appeared (see Fig. 1 below). For 
small lattice dimensions, they observed the convergence of this phase to a fixed 
point and concluded that asynchrony had a stabilizing effect on Life [3]. 

The phase transition was then measured with precision by Blok and Bcrg- 
ersen, who used the final density (i.e, the fraction of l's sites) as a means of 
quantifying the phase transition. They measured the value a c for which the 
phase transition was to be observed and found a c = 0.91 [5]. They showed that 
the type of phase transition is continuous (or a second-order transition): when a 
is decreased from a = 1.0 to a = a c , no change is observed in terms of the values 
of the average density. When we have a < a c the "labyrinth phase" gradually 
appears and the average density starts increasing in a continuous way. It is thus 
the derivative of the density that shows discontinuity rather than the function 
itself. 

Figure 2 shows that the removal of links between cells does not qualitatively 
perturb the aspect of the final configurations attained. So, according to the 
observation of steady states, synchronous Life seems somehow robust to topology 
perturbations. However, we also noticed that the transients are much shorter in 
presence of topology errors: for N = 50 x 50, the order of magnitude of transients 
are T = 1000 for e~ = and T = 100 for e" = 0.1. 




a = 1.0 a = 0.5 sequential updating 



Fig. 1. Life configurations for N — 50 x 50, after T — 100 time steps, starting from 
a random configuration of density di n i = 0.5. In the sequential updating, cells are 
randomly updated one after another. 
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Fig. 2. Life configurations for synchronous evolution (a — 1.00) with N = 50 x 50, 
after T = 1000 time steps, starting from a random configuration of density dim = 0.5. 




<T = e - = 0.05 e" = 0.10 



Fig. 3. Li/e configurations for TV = 50 x 50, after T = 1000 time steps, starting from a 
random configuration of density di n i = 0.5: (up) a = 0.90 (middle) a = 0.75 (bottom) 
a = 0.50. The figure in the upper-left corner shows that the system is still in a transient 
mode. 



Figure 3 shows what happens when both asynchronism and topology pertur- 
bations are added. Rows of Fig. 3 display the behavior with a fixed synchrony 
rate a and columns display the behavior with a fixed missing- link rate e~ . We 
see that increasing topology errors from e~ = to e~ = 0.05 makes the phase 
transition occur for a higher value of synchrony rate a c . With a further increase 
from e~ = 0.05 to e~ = 0.10, the phase transition cannot be observed any more, 
at least for the selected values of a. 

This demonstrates that both parameters e~ and a control the phase tran- 
sition between the "inactive-sparse phase" [9] and the "labyrinth phase" (LP). 
The next section is devoted to quantitatively measure the interplay of these two 
control parameters. 

3.2 Quantitative Approach 

To detect the apparition of the labyrinth phase (LP), we need to look at the 
configurations by eye or to choose an appropriate macroscopic measure. Clearly, a 
configuration in LP contains much more l's than a configuration in the "inactive- 
sparse phase" ([9]). This leads us to quantify the change of behaviour using 
the measurement of the "steady-state density" (i.e. the average density after a 
transient time has elapsed). This method has been chosen by various authors 
(e.g. Blok and Bergersen [5]) and it has been applied to exhaustively study both 
the dynamics [7] and the robustness to asynchronism [8] of one dimensional 
elementary cellular automata. 

We define the steady-state density p(di n i,a) using the sampling algorithm 
defined in [8]: Starting from a random initial configuration constructed with 
a Bernoulli process of parameter <ii n i, we let the AC A evolve with a synchrony 
rate a during a transient time Transient then we measure the value of the density 
during a sampling time T samp iing- The value of p is the average of the sampled 
densities. 

The sampling operation results in the definition of a function p(d; n i, a) that 
can be represented in the form of a "sampling surface" . This surface contains part 
of the information on how the behaviour of a CA is affected by asynchronism. 
Figure 4 shows the experimental results obtained for N = 50 x 50, T tran siont = 
1000 , T samp ii ng = 1000. The transient time Ttransiont is chosen according to the 
observations made in [1] where equivalent transient times were found for greater 
lattice sizes. 

Let us first look at what happens for d; n i £ [0.2, 0.8] (right column of Fig. 4). 
The invariance of the surface relatively to the d; n i-axis shows that the macro- 
scopic behavior of Life does not depend on the value of this parameter within 
this range. The upper right corner of Fig. 4 shows that for e~ = (regular 
topology), the phase transition occurs for a c ~ 0.90 as expected [5]. However, 
when e~ increases, experiments show that a c also decreases. This means that the 
settlement of LP becomes more difficult as links are removed; this can be inter- 
preted as an increase of the robustness to asynchrony. We can observe that for 
e~ = 0.10 , the surface is flat and horizontal, which means that the behavior is 




Fig. 4. Sampling surfaces for e~ = [0---0.10], N = 50 x 50, Transient = 1000 , 
Trampling = 1000. The left column has a different range for d ln i to focus on the be- 
havior for small initial densities. 



not anymore perturbed by asynchronism (at least if we consider our observation 
function). 

The left column of Fig. 4 shows the behavior of Life for di n ; £ [0,0.2]. We 
observe two different abrupt change of behaviors. On the one hand, there is a 
value of df ni which separates the "inactive-sparse phase" and LP. On the other 
hand, the value of df n[ increases as e~ increases. This means that LP becomes 
more difficult to reach when links are removed; which again can be interpreted 
as a gain of robustness. 

Experiments were held for various lattice sizes and allowed to control that 
the sampling surface aspect was stable with TV; we however observed that when 
e~ and a are fixed, the value of <i? ni is a decreasing function of N. 

All the previous phenomena may be the consequence of multiple intricate fac- 
tors. In the next section, we study the evolution of so called "micro-configurations" 
put in an empty array and propose a first hypothesis in the direction of under- 
standing these behaviors. 

4 Micro-configurations Analysis 
4.1 Experiments 

The observation of the settlement of LP shows that it can develop from very 
localized parts of the lattice and then spreads across the lattice until it fills it 
totally. By analogy with a crystal formation, we call "germs" these particular 
configurations that have the possibility to give birth to an invasive LP. We inves- 
tigate the existence of germs by performing an exhaustive study of the potential 
evolution of micro- configurations, i.e. 3x3 configurations that are placed in an 
empty array. There are 512 such configurations and we experimentally quantify, 
for each one, the probability that a it becomes a germ. Our goal is to infer the 
behavior of the whole structure from the evolution of these micro-configurations. 

Setting the synchrony rate to a = 0.5, we used the following algorithm: 
For every micro-configuration i £ I, (a) we initialize the lattice N x N with i, 
and (b) we let the CA evolve until it reaches a fixed point or until it reaches 
LP. We repeat S = 1000 times operations (a) and (b) for the same initial micro- 
configuration but for a different update histories. We consider that the CA has 
reached LP if the density is greater of equal than a limit density c?oo = 0.1. 
Indeed, we observed that if the CA was able to multiply the number of l's from 
the micro-configuration to a constant ratio of the lattice, then it will almost 
surely continue to invade the whole lattice and, asymptotically, reach LP. We 
experimentally obtain the probability -P ge rm[i] that a configuration i is a germ. 
Grouping micro-configurations by the number k of l's they contain, we obtain 
an array with 9 entries P ge rm(k), k £ [0, 9], displayed in Table 1. 

Results show that for k < 3, P gC im(k) = 0, which means that all such micro- 
configurations always tend to extinction. For n > 3, the probability to reach LP 
increases as n increases. 



Table 1. Probability P gcrln (k) for a micro-configuration to be a germ (in %) as a 
function of the number k of l's for different missing- link rates e~ (in %). 
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Fig. 5. Probability to reach LP as a function of di n i for different lattice sizes N £ 
{20, 30, 40, 50} and different missing-link rates: <T — 0, <T = 0.02, <T = 0.04 from top 
to bottom. Left column shows theoretical curves calculated with the independent-germ 
hypothesis; the right column shows the actual measurements. 



4.2 Inferring Some Aspects of the Global Behavior 

Our idea is to use the previous results about germs to extend them to a descrip- 
tion of the global behavior of Life. Unfortunately, we are not able to do that 
in the exact way. In a first approximation, let us assume that we can infer this 
global behavior by approximating the probability Plp for a uniformly initial- 
ized system to reach LP using an "independent-germ hypothesis": interactions 
between potential germs are neglected and we assume that LP is reached if and 
only if there is at least one cell that gives birth to LP. This results in the appli- 
cation of formula : Plp = 1 — (1 — Plpi) N * N where Plpi is the probability that 
one cell gives birth to LP. We have Plpi = X)fe=o Pgerm(k) 1 Papp(fc, d- ln i), where 
Pa PP (fc,rfini) is the probability that a micro-configuration initialized randomly 
with di n i contains k l's. It is simply obtained by applying the binomial formula: 

P (k rl- ■) — ( 9 )rl- - k (\ - (]■ -) 9 ~ k 

Calculated and experimental values of P\jp{d- m \,e~) are given in Fig. 5 and 
show that this assumption is justified in a first approximation even if the pre- 
dictions seem more accurate for small values of e~ . 

The germ hypothesis allow us to understand better some of the observed 
behavior. Let us first consider the abrupt change of behavior observed for di ni e 
[0,0.20]: this can come from the fact that as <ii n i increases the probability to 
observe a micro-configuration that contains more l's increases thus increasing 
the probability to find a germ in the initial configuration. We can also understand 
with this point of view the invariance of the sampling surfaces in the <i; n i-axis 
with d[ n [ > 0.20 by the fact, observed in Fig. 5 that in this case Plp ~ 1, that is 
the "labyrinth phase" always appears. In the same way, the shift of df ni observed 
in Fig. 4 when varying e~ can be explained by the looking at the variations of 
Pgerm(fc) with e~ : we see that all probabilities to reach LP decrease when e~ 
increases. Finally, we are able to qualitatively predict the scaling of df ni with 
the lattice size N from the plots the function PLp(<iini, e - ) : when <i; n i and Plp 
are relatively small (i.e., Plp < 0.1), we have a linear scaling of Plp with N 2 ; 
whereas as Plp is close to saturation, there tends to be no variations with N. 

5 Conclusion 

Experiments have shown that Life's transition from an "inactive-sparse phase" 
to a "labyrinth phase" (LP) is a continuous phase transition dependent on the 
synchrony rate a and whose critical value a c is controlled by the missing-link 
rate e~ . As the topology was perturbed (i.e., when e~ increased), the inactive- 
sparse phase domain extends while the LP domain shrinks. The abrupt change 
in behavior according to the values of initial density dj n j was interpreted with the 
hypothesis that the settlement of LP results from the development of "germs" , 
i.e. small configurations that are able to "colonize" the whole lattice. The study 
of the evolution of potential germs from micro-configurations allowed us to start 
to understand the observations and to give some predictions on the probability 
to reach LP starting from a random configuration. One interesting question is 



to examine whether these observations hold for a large class of CA or if they are 
somehow related to the computational universality of Life. 
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